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ABSTRACT 

The HDO/H 2 abundance ratio is thought to be a key diagnostic on the 
evolution of water during the star- and planet-formation process and thus on 
its origin on Earth. We here present millimeter-wavelength high angular reso- 
lution observations of the deeply embedded protostar NGC 1333-IRAS4B from 
the Submillimeter Array targeting the 3i2 — 221 transition of HDO at 225.6 GHz 
(E u = 170 K). We do not (or only very tentatively) detect the HDO line toward 
the central protostar, contrasting the previous prominent detection of a line from 
another water isotopologue, H^O, with similar excitation properties using the 
IRAM Plateau de Bure Interferometer. The non-detection of the HDO line pro- 
vides a direct, model independent, upper limit to the HDO/H 2 abundance ratio 
of 6xl0 -4 (3a) in the warm gas associated with the central protostar. This upper 
limit suggests that the HDO/H 2 abundance ratio is not significantly enhanced 
in the inner w 50 AU around the protostar relative to what is seen in comets and 
Earth's oceans and does not support previous suggestions of a generally enhanced 
HDO/H 2 ratio in these systems. 



Subject headings: astrochemistry — stars: formation — protoplanetary disks - 
ISM: abundances — ISM: individual (NGC 1333-IRAS4B) 
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Introduction 



One of the key questions concerning the formation of planets, and our own Earth in 
particular, is how water has evolved on its way from collapsing molecular cloud cores, to 
protoplanetary disks and eventually Earth's oceans. In c old and quiescent reg i ons, the gas- 



phase water abundance is low, only 10~ 9 — or less (e.g.. lBergin fc Snellll2002l : ICaselli et al. 
20101 ). but in regions forming low-mass protostars with intense heating or active shocks, its 
abundance can reach 10 -4 with respect to H2 - comparable to or higher than that of CO 



(e.g., iKristensen et al.ll2010l ). Thus, in the earliest prestellar stages and in large parts of the 
dense envelopes around low-mass protostars, H 2 is largely frozen out and the dominant 
constituent of the icy mantles of dust grains. It is thought that the amount of deuterated 
relative to non-deuterated water is established in these grain mant les. Observations show 
that this abundance ratio is enhanced in comets (~ 3 — 4 x 10~ 4 ; e 

1.5 x 10~ 4 



Bockelee-Morvan et al. 



Robert et al 



2000 



and 



19981 ; IVillanueva et all 120091 ) and Earth's oceans (~ 1.0 x ill"; e.g. 
references therein) relative to primor dial D/H ratio in the protosolar nebula (~ 1.5xl0 -5 ; 
Geiss fc Gloecklerlll998l ; lLinskyll2003l ). Determining when and where this ratio is established 
is thus a potentially important way to trace the evolution of water in the different stages of 
star formation. In particular, it might help determine what fraction of water was brought 
to Earth at later times through c ometary impacts compared to being accrete d during the 
formation of the planet itself (e.g.. iMorbidelli et al.ll2000l ; iRaymond et al.ll2004l ) and whether 
water undergoes signficant processing in warm regions of protostellar envelopes or disks 
before reaching newly formed planets. 

Previous attempts at measuring the ratio of HDO and H 2 in low-mass protostars 
have had mixed success. Direct attempts to detect HDO in the solid phase at infrared 
wavelengths have pro vided upper limit s of [HDO]/[H 2 0] varying from 0.005 to 0.02 in a 
sample of protostars (IParise et al.l 120031 ). which are not strong constraints compared to the 
two orders of magnitude lower abundance ratios observed in, e.g., the cometary systems 
about two orders of magnitude lower. Measurements at (sub)millimeter wavelengths could 
potentially be sensitive to lower abundances of HDO in the gas-phase in regions where the 
ices have evaporated off the dust grains and thus trace their compositions indirectly, but 
such measurements are complicated due to the achievable angular resol ution using single- 
dish radio telescopes an d have led to discrepant results. As an example, IStark et al.l (j2004f ) 
and IParise et al.l (120051 ) analyzed the HDO emission toward the deeply embedded protostar, 
IRAS 16293-2422 compared to previous ISO and SWAS observations of non-deuterated H 2 0. 
Both pr esented detai led models for the HDO abundance structure with significantly different 
results: IParise et al.l found evidence for a drop in the [HDO]/[H 2 0] ratio from 0.03 in the 
hot gas in the inner envelope (i.e., significantly e nhanced abo ve the cometary values) to 
<0.002 in the cold gas on larger scales, whereas IStark et al.l in contrast found that the 
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abundance of HDO was best described as being constant throughout the envelope and with 
a [HDO] / [H 2 0] ratio of 2 x 10~ 4 , comparable to the cometary value, in the part of the envelope 
with T > 14 K. Part of this discrepancy may arise due to the complexities in comparisons 
between ground-based single-dish observations of HDO and space-based observations of H 2 
probing signicantly different spatial scales - and lines excited in different regions of the 
protostellar environment. 

Our detection of thermal H^O emission toward the deeply embedded Class low-mass 
protostar, NGC 1333-IRAS4B (d = 250 pc), with the IRAM Plateau de Bure Interferometer 
(jJorgensen fc van Dishoeckl 120101 ) makes it possible to revisit the discussion of the origin of 
water on Earth - namely where the relative abundances of deuterated and non-deuterated 
are established. In this letter we present a search for the 225.9 GHz 3i 2 — 2 21 transition of the 
HDO isotopologue at high-angular resolution (2.5-3"; corresponding to 625-700 AU) using 
the Submillime ter Array (SMA). This transition is close in energy to the H 2 8 3 



13 
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transition from I Jprgensen fc van Dishoeckl (120101 ) (E up = 170 K vs. 204 K) and a straight- 
forward comparison can thus be made between these species with minimal uncertainties due 
to their excitation properties. 



2. Observations 



W e observed NGC 1333-I RAS4B (IRAS4B in the following; a=0 3 h 29 m 12.0O: 1 £= +3i d 13'08"l 
[J2000] IjQrgensen et al.ll2007h using the Submillimeter Array (SMA; IHo et alJliopj ) on 2008 
October 31. The receivers were tuned to the HDO 3i 2 — 2 21 transition at 225.89672 GHz and 
the SMA correlator assigned a chunk with 512 channels centered on this line corresponding to 
a 0.27 km s _1 spectral resolution over about 100 MHz of the chunk. The chunk was placed 
so that the upper sideband covered the D 2 CO transition at 236.102128 GHz with similar 
spectral resolution. For the remaining SMA passband 128 channels were allocated in each 
chunk giving a uniform spectral coverage of 1.1 km s -1 from 225.0 GHz to 227.0 GHz (LSB) 
and 235.0 to 237.0 GHz (USB). The source was observed in the compact configuration of 
the SMA, which at the time of observations had 7 antennae operational providing baselines 
ranging from 7 to 60 kA. 

The data were calibrated using the SMA version of the MIR package and imaged using 
Miriad. The calibration followed the standard approach: the absolute flux calibration was 
established through observations of Uranus, the bandpass by observations of the strong 
quasar 3c454.3 and the complex gains by regular observations of the nearby strong quasars 
3c84 and 3clll (approximately 5.2 and 4.7 Jy at 1.33 mm, respectively). Integrations with 
clearly deviating amplitudes and/or phases were flagged and the continuum was subtracted 
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prior to Fourier transformation of the line data. With natural weighting the resulting beam 
size is 3.1"x 2.5"at a position angle of -66°; the field of view is 52" (HPBW) at 1.33 mm. The 
resulting RMS noise level is 58 mJy beam -1 channel -1 for the line data with the 0.27 km s -1 
spectral resolution using natural weighting. 



Results 



Both IRAS4B and its nearby companion IRAS4B' are clearly detected in the continuum 
data. Table [T] lists the results of elliptical Gauss ian fits to the two sourc es: both are resolved 
with fluxes in agreement with the results from Ijprgensen et al.l (120071 ) within the absolute 
flux calibration uncertainty of 20%. 

Fig. [1] shows the spectrum from 1 GHz of the LSB of IRAS4B labeling a number of 
lines with peaks at the 3a level or above. The H 2 CO 3i 2 — 2 1;1 transition at 225.698 GHz 
is clearly dominant as expected, with lines of complex organic molecules as well as OC 34 S 
and C-C3H2 also detected. Most importantly the identification of the multiple lines centered 
at the continuum position confirms both the pointing of the observations as well as the 
calibration of the frequency axes of the spectra. 

Fig. [2] shows a zoom-in on the spectrum in the high resolution chunk around the HDO 

31.2 — 22,i transition. No clear line is seen at the system ic velocity of 7.0 km s -1 of the H^O 

31.3 — 2 2 transition ( Ijprgensen fc van Dishoeck! 120101) a nd other species surveyed at high 
angular resolution with the SMA ( Jprgensen et al.l2007 ) - although a tentative detection 
may be present at the ~3cr level (integrated intensity). This is confirmed by the integrated 
intensity map (Fig. [21 right), which shows a 3a contour close to the location of the continuum 
peak. Statistically speaking it cannot be counted a firm detection, however, and we therefore 
treat it as an upper limit in the following. 



Table 1: Parameters for IRAS4B and IRAS4B' from elliptical Gaussian fits to their continuum 
emission. 





IRAS4B 


IRAS4B' 


Flux 


1.021 Jy 


0.364 Jy 


R.A. (J2000) 


03:29:12.0 


03:29:12.8 


DEC (J2000) 


31:13:08.1 


31:13:06.9 


Extent 


1.7"x0.96" (-68°) 


1.6"x0.91" (+55°) 



a Size of Gaussian from fit in (w,t>)-plane (i.e., deconvolved FWHM size) and position angle 

of major axis (in parentheses). 
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4. Discussion: an upper limit on the ratio of the HDO/H 2 abundances in 
the inner 50 AU regions of a low-mass protostar 

The 3cr upper limit on the integrated line intensity (or approximate strength of the 
tentative detection) of the HDO transition can be estimated as 3o"i nt = 3v Av 5v a Tms « 
0.13 Jy km s _1 , assuming a line width to zero intensity of 2 km s _1 as found for the H^O 
transition. The column density of HDO corresponding to this intensity can be estimated 
under the same assumptions as for the H^O transition - and can thus be compared directly, 
assuming that the emission is optically thin and uniform over the extent of the H^O emission 
(thus diluted in our SMA beam). 



Using the non-LTE escape probability code Radex (Ivan der Tak et al.l 120071 ) we calcu- 
lated column densities corresponding to the observed line intensities for different values of 
the kinetic temperature and H_2 density. Fig.[3]shows the derived HDO/H2O abundance ratio 
as function of kinetic temperature and density, respectively. In the range from about 100 
to 1000 K and for densities larger than about 10 9 cm -3 (the density expected in the inner 
envelope or disk) the derived [HDO]/[H20] 3a upper limit on the abundance ratio varies 
from about 2xl0~ 4 to 6xl0 -4 assuming an ortho-para ratio for H^O of 3:1 and a 16 0/ 18 
ratio of 560. If the tentative detection would be confirmed by more sensitive observations, 
the implied HDO/H 2 abundance ratio would suggest a level of deuteration comparable to 
what is observed in cometary ices. The Radex calculations show that the inferred abundance 
ratio can vary from about lxl0~ 4 to lxl0~ 3 if the excitation temperatures are allowed to 
differ for the two species and vary between 100 and 1000 K, covering a wider range of phys- 
ical parameters. The beauty of the above result is that it does not depend on the assumed 
physical structure of the source - and in particular, whether the origin of the observed water 
is due to shocks or ice desorption in either envelope or disk material. 

One caveat is that the H^O transition could b e masing as is the case for the same 



transition of main Hj, 6 isotologue at 183 GHz (e.g., ICernicharo et al.lll994l ). In that sce- 
nario the H^O column density would be overestimated and thus the [HDO]/[H 2 0] ratio 
underestimated. A few observational facts argue against this being a major issue for the 
H^O transition toward IRAS4B, though: the width of the Hj, 8 transition (2 km s~ x to 
zero-intensity) is very close to those of the other species (SO2, C2H5CN and CH3OCH3) 
observed in the central beam toward the continuum peak of IRAS4B. This trend is further 
strengthened with detections of the Hj, 8 transition toward two other embedded protostars 
(M. Persson et al., in prep.): the lines of Hj, 8 0, SO2 and the organic species in these three 
sources show smaller species-to-species variations than source-to-source variations both in 
terms of their line widths and strengths. This suggests that the origin of the emission is the 
same in each of the individual sources and therefore that H^O emission is dominated by 
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the thermal component and at most weakly masing. Finally, previously detect ed H2O water 



masers at 22 GHz in IRA S4B are clearly offset from the systemic velocity (IMarvel et al. 



20081 ; iDesmurs et al.l 120091 ) and their motions in the outflow propagation direction are per- 



pendicular to the tentative velocity gradient observed in the H^O transition. Thus, the 
H^O line is at least not masing in the gas where the 22 GHz H2O maser has its origin. 

The d erived upper lim i t to t he H DO/H?0 abundance ratio does not support the sug- 
gestion by ICeccarelli et al.l (120051 ) and iParise et al.l (120051 ) that the HDO/H 2 abundance 
ratio is abnormally high in low-mass protostars or their disks compared t o the abundance 



ratios found in comets, Earth's oceans or in high -mass hot core region s fe.g..|Jaca et al 



Gensheimer et al.lll996l ; Ivan der Tak et al.ll2006l ). In agreement with lStark et al.l (12004 ). the 



1990 



results favor a scenario in which the HDO/H 2 ratio is conserved during the star forma- 
tion process. There is thus also no need to invoke the formation of water in the gas-phase 
in a longer-lived early stage of low-mass protostars where the deuterium fractionation is 
high to explain the constraints on the HDO/H2O abundance ratio. The conclusion that the 
HDO/H 2 abudance ratio is determined in the cold core phase and is conserved throughout 
the formation of l ow-mass star s and icy solar-system bodies is consistent with the evolu- 
tionary models of IVisser et al.l (120091 ) that show that the bulk of the ice remains in solid 
form during infall and disk formation. It is possible that the physical mechanisms for the 
origin of the warm gas in NGC 1333-IRAS4B and IRAS 16293-2422 differ: the latter source 
shows a complex structure with multip le organic species peaking at different locations (e.g., 
Kuan et al.ll2004l; iBisschop et al.ll2008f) and likely interaction with the protostellar outflow on 



small scales (I Chandler et al 



20051 ). SMA observations of that source (J. K. J0rgensen et al., 



in prep.) suggests that this complex interplay could also have an impact on the line emission 
for the HDO isotopologue underlining the need for high angular resolution observations. 

In any case, these observations have shown the great potential offered by high angular 
resolution submillimeter (interferometric) observations, even from the ground, for tracing 
the origin and evolution of water vapor during the evolution of young stars. An important 
step forward would be to spatially resolve both the emission profiles of HDO and H^O in 
protostars in different evolutionary stages to reveal any temporal or spatial variations of their 
abundance ratios. More sensitive observations with the Plateau de Bure as well as future 
observations with the Atacama Large Millimeter Array in Band 5 will make this possible and 
thereby make an important contribution toward answering the questions about the origin of 
water in planetary systems. 



We are grateful to Tim van Kempen for helpful comments about the manuscript. Re- 
search at Centre for Star and Planet Formation is funded by the Danish National Research 
Foundation and the University of Copenhagen's programme of excellence. Research in as- 
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Fig. 1. — Spectrum extracted in the central 3.0" x 2.5" beam toward the continuum position 
for IRAS4B. The detected lines are indicated at the position of their catalog rest frequency 
corrected for the 7.0 km s _1 systemic velocity of IRAS4B. The box shows the zoom-in around 
the HDO 3i,2 — 2 2 i transition displayed in Fig. [21 the high spectral resolution chunk including 
this line has been rebinned to the same 1.1 km s -1 spectral resolution as for the other chunks 
in the wider bands. The horizontal dotted lines indicate the zero- and 3a levels. 
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Fig. 2. — Left: Spectrum centered on the expected location of the HDO 3i )2 — 2 2) i transition. 
The dotted lines indicate the zero- and la levels. The dashed line, the expected location of 
the HDO line at the systemic velocity of the H^O 3^3 — 22,o transition and other species at 
7.0 km s _1 . Right: Integrated intensity of the data cube over ± 1 km s _1 around the expected 
location of the HDO transition. The contour corresponds to the 3a level (no emission is seen 
at 4er levels or higher). 
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Fig. 3. — Estimates of the HD0/H 2 abundance ratio as function of temperature (left) and 
density (right) from non-LTE Radex calculations. Each panel shows three curves: in the left 
panel for constant densities of 10 9 cm -3 , 10 10 cm~ 3 and 10 11 cm~ 3 and in the right panel for 
constant temperatures of 170 K, 300 K and 600 K. 



